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ABSTRACT 

We evaluate the inclusive semileptonic B decay spectrum as a function of the final-state 
hadrons energy to second order in the inverse quark mass expansion and tree level in as- We 
argue that there is an energy interval below the c production threshold that could be used to 
determine Vub- 



March 1994 



^Internet: bouzas@physics.ucla.edu 
^ Internet : zappala@physics.ucla. edu 



1. Introduction. The inclusive semileptonic decay of the B meson is of great importance for 
the determination of the Cabibbo-Kobayash-Maskawa (CKM) matrix elements Vcb and Vub- For 
this reason several attempts have been made to predict the total width and the electron energy 
distribution in this decay. It has recently been shown |j^ that both can be computed in the 
framework of the QCD operator-product expansion, without resorting to models for the structure 
of the initial hadron. To lowest order the parton model result is recovered. Matrix elements of 
higher dimension operators are then expanded in powers of the inverse of the b quark mass 1/m by 
applying heavy-quark effective theory (HQET) Q. 

This method has been applied |^ ^ (see also ^) to compute the electron energy spectrum to 
order 1/m?. However, for values of the electron energy in the vicinity of its maximum kinematical 
limit the corrections to the parton model are unreliable, being dominated by large unphysical 
oscillations and furthermore singular at the endpoint. This fact spoils the possibility of determining 
Vub^ since this quantity must be extracted from the endpoint region where the b ^ c channel is 
kinematically forbidden. Some tentative solutions to this problem have been suggested such as 
gaussian smearing or integration of the spectrum in the endpoint region. Also a form factor 
formalism has been proposed Q. 

In this letter we compute the inclusive semileptonic differential decay width as a function of 
the total energy of the final-state hadrons employing the approach outlined above. We suggest 
that this quantity could be suitable for the determination of Vub- In this case, the relevant part 
of the spectrum corresponds to total hadronic energies smaller than the charm mass where b ^ u 
is the only decay mode available and it is worthwhile to study the behaviour of the heavy-quark 
approximation scheme within this region. 



2. Hadronic spectrum. In order to obtain the hadronic energy distribution we closely follow 
the procedure used in |^ for the electron spectrum. For that reason, we omit all technical details 
related to the computation. In the sequel we neglect radiative corrections of order Qs/tt. 

We start from the weak hamiltonian 



n{x) = -v,b^Jt{x)Mx) (1) 



where Gp is the Fermi constant, Vjb is the CKM matrix element for b ^ j and Jj^ and are the 
leptonic and hadronic weak currents respectively, 

Jtix) = Mx)r{l-l5)Mx) (2) 



Jti{x) = ^/'j(x)7^(l - 75)V'6(x) (3) 

The squared amplitude for the semileptonic B decay can be written as {k and k' are the electron 
and neutrino momenta respectively, and q = k + k'), 

\A\^ = \3]^^2n)L>^^{k,k')W,A<l) (4) 
The leptonic tensor L^^, appearing in ^ is given at tree level by the expression]^, 
$^(/,F|jr(x)|0)(0| J/'^(0)|/,F) = e"^-L'^''ik,k') 

pol 

= e^''^ 8 i^k'^k"' + k'^k'' - {k ■ k')g^"' - ie^"'''^kak'p) (5) 

The hadronic tensor, defined as, 

Wf^'^il) = h^J d'xe^''{B\4mX){X\Mx)\B) (6) 

(where the sum runs over a complete set of hadronic states) can be related to the discontinuity of 
the matrix element of the time ordered product of hadronic weak currents, 

T^u = -ij d^xe'^^B\TJl{Q)Ux)\B) (7) 

through the relation, 

Im(VL^^) = -7rT^^,L'^^ (8) 
Finally, using (^) and (^), the B decay width is given by 

Tb = \Vj,\'Gl J ^^^6{k')6{ko)6{k")eik'o)lm{-T,,Ln (9) 

with the normalization {B{Mbv)\B{Mbv')) = {2TTfv^5'^{MBV - Mbv') [§, Mb and V being the 
mass and velocity respectively of the B meson. The problem of calculating the total decay width 
and the leptonic and hadronic spectrum is then related to the evaluation of the tensor T^*^. We 
shall limit ourselves to the case j = u and the u quark mass will be neglected. 

Due to the T-ordering the matrix element in contains a u quark propagator. In the following 
we shall treat separately two contributions to the quantity T^^i^L'^'^ 

T^^Lf'" = tWl^'^ + r^2)^A«- (10) 



^We use the convention £9123 = 1- 
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where the first term contains the lowest order approximation for the u propagator [i.e. the free 
propagator) whereas in the second one the first order correction to the propagator, consisting of a 
gluon insertion, is taken into account. 

The first step to compute tI:}J L^y is to factorize out the large x-dependence in the b field, 

V^6(x) = e-^"^^^V'6(x) (11) 

with m being the b quark mass and v the 4-velocity of the B meson. The field on the right-hand side 
of (|ll]) can be split into two components ^^(x) = h{x)+x{x) defined by the properties ^/i(x) = h{x) 
and '^x{x) = ~x{x)- xi^) can be eliminated in terms of h{x) through an expansion in powers of 
1/m 0. 

At the same time, we can perform a short-distance expansion since the x dependence of ip'i,{x) 
is small in the sense that it is related to the "residual" momentum k, defined as, 

K = pb — mv (12) 

K is non zero for finite m because the b quark does not carry all of the B momentum, but it vanishes 
in the limit of very large quark mass. 

The combination of short-distance expansion and HQET leads to a 1/m expansion of the matrix 
element in terms of operators of increasing dimensions. The lowest order coincides with the the 
semileptonic decay of a free b quark, since it corresponds to setting k = 0, (i.e., pb = mv) and thus 
identifying the decaying particle with the b quark. Corrections of order 1/m vanish ||l|, whereas 
those of order 1/m^ are expressed in terms of the two matrix elements ^ 

Kb = -{Bm^-^hmB) (13) 

Gb = Zb{B\-hiO)^^^^h{0)\B) (14) 

D being the QCD covariant derivative, g the strong coupling constant, ctq,^ = (i/2)[7Q,, 7^], and 
ZbG^^ the renormalized gluon field strength tensor. The quantity Gb is directly related to the 
experimentally measured mass splitting M{B*) — M{B) from which one gets the small value Gb = 
—0.0065. The numerical value of Kb is estimated to be Kb = 0.01 ^ from QCD sum rules (see 

i)- 

The explicit expression for T^h) L^^^ to the order l/m^ is given by, 

r^L^- = 641 d^x J ^e-*(-''-'?-f)--2^(i?o + Ri + R2) (15) 



where p is the n-quark momentum and, 

Ro = {p- k) {v • k') 



R2 



The discontinuity of T^ll^L^"' in (|T|) is given by the u quark propagator 

1 



^{vk') vpyp-k) 
^m'^Kbip -k) {v k') (x^ - {v ■ xf) 



(16) 



Im 



+ ie 



-tt6{p') 



The integration over x in (15) yields terms proportional to 5'^{mv—q—p) and its first two derivatives, 
according to the power of x appearing in Ri. 

We notice that each term Ri of the expansion is symmetric in the two momenta p and k. This 
symmetry can be related to a Fiertz identity for the product of currents in the amplitude that 
exchanges the electron and quark fields ipe-, i^u- Thus, from (^), (|T5|) and (|l6|), the first piece of 
the decay width is invariant under exchange of p and k and the two spectra dT^^ /dEf. and 
dV^^ /dEu have exactly the same form. 

After performing the integrations we obtain. 



1 dV 



(1) 

B 



F dw 



20 n n 



e{l-w) 



+ 



5{l -w) + -Kb5'(l 



w 



(17) 



where we introduced the dimensionless variable w = 2Eu/m and the lowest order width F = 
(|KfepG2,m5)/(192^3), 

Let us now consider the second term TjiJh^^ in (0). The complete expression for T^^^L'^'^ to 
order 1/m^ includes the correction obtained by attaching a gluon line to the free u propagator. 
The insertion of the gluon field on the propagator in the matrix element of (0) generates a term 
proportional to Gb Q . Notice that this insertion is not related to the one loop self-energy correction 
to the propagator, whose imaginary part is proportional to the probability of real gluon emission. 
The new term is, then. 



d^x 



d'^p 



-i(mv~q- 



~p)x(_}_\ 

Vp2 -\- is) 



m^GbX 



[{v -k) {p- k') - {k ■ k') {v ■ p)] 



(18) 



Unlike Tjh) L^^, in the above equation the expression in square brackets is antisymmetric in p and 
k. 
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In fact the presence of the gluon field in the hadronic current breaks the symmetry between 
ipe and ipu- In this case one can stih apply a Fiertz transformation exchanging -i/^e and the string 
ipu{ipuC^i^u)- This operation yields again the result ([l8|), but it does not simply correspond to the 
exchange of p and k in the scalar products of (pisD- It also has the effect of substituting —Gh for 
Gb- 

In (|l8|) the discontinuity generated by the squared propagator is, 

Im(^^)'=7r5'(p2) 

Due to the integration over po, the derivative acting on the delta function can be transferred to 
the other terms in (^). This way lm{T^ L'^'^) contains no derivatives of 5{p'^) but the derivative 
of the momentum conservation delta function has been generated. Thus, we can see that both in 
the first part of the calculation in (^) and in the gluon emission part ([Tsl), the potentially singular 
terms appear at this step of the computation as derivatives of 6'^{mv — q — p). 

The contribution of T^u L^'^ to the spectrum, obtained by integrating in (|18|) is. 



1 r/r*-^^ 4 r 1 

^-j^ = -Gb\{4w - 3w^)9{l -w)- 6{1 - w)] (19) 
1 aw 6 L J 

As already mentioned, the gluon contribution to the leptonic spectrum has a different form. From 
dH) one gets {y = 2Ee/m), 



r dy 3 



Gb -9y' 9(1 -y) + 35(1 - y) (20) 



Incidentally, we notice that both (19) and (P0|), although relevant in the decay energy spectra, give 
zero contribution to the total B decay width. 

Finally, the u quark decay spectrum to the order 1/m^ is given by the sum of ( [l7| ) and (19) 
and it is plotted in Fig. 1. Due to the presence of the delta functions we have limited the range of 
w to < w < 0.9 in the figure, and to show the small Kb dependence curves for Kb = 0, and 0.05 
(dotted curves) are also included. At w = 0.5 we find that the differential width changes by about 
5% as Kb varies from to 0.05. 



3. Discussion. The differential decay rate for the process B eVX is obtained, as shown above, 
from the discontinuity of T^^L^i, whose cut is given by the propagator of the final u quark (or c 
quark if X contains charmed mesons) only. It is natural, then, to identify the total energy of the 
final-state hadrons Ex with the energy carried by the u (c) which is the only particle contributing 
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to the imaginary part of the ampUtude. This would obviously not be true had we included 0(as/vr) 
radiative corrections. In this case the propagator would be dressed by self-energy corrections, so 
that for instance a real gluon emission would give a new contribution to the discontinuity. 

Let us now briefly comment on our results. Equations ( [T7| ) and (|l9|) show that the hadronic 
spectrum has the same kind of singularities at the endpoint w = 1 that characterize the electron 
spectum at y = 1. To explain this point one should notice that both hadron and electron spectra 
for the decay of a free b quark (that is our lowest order approximation) do not vanish smoothly at 
the endpoint w = 1 and y = 1 but contain functions 6(1 — w) and 9{1 — y). Singularities appear 
in the following corrections due to expansion of those step functions ||3|, ^] . More explicitly, looking 
at (^) in momentum space, we just have the expansion of 6^{pb — p — q) around k = 

5^{pi, — p — q) = S^{mv — p — q) + n ■ d5^{mv — p — q) + ... (21) 

and performing the integration, due to the symmetry of the variables p and k, the derivatives in ( ^l]) 
are transformed in derivatives 6{1 — w) when computing cLTb/ dw or of 5(1 — y) for dV b /dy, thus 
generating analogous singularities for the two cases. The same argument is valid for the singularity 
in the gluon insertion correction (|l8|) since, as already discussed, it can be expressed again as a 
derivative on 5'^{mv — p — q). 

Another point to be discussed concerns the region around w = 0. As explained in IQ the 
discontinuity of T^'^L^j^y relevant to the problem considered here is given, in the {v ■ q) complex 
plane at fixed q"^, by a cut along the real axis in the range < v ■ q < (M^ + q"^ — MI)/{2Mb) 
(M/j is the mass of the lightest hadron that can be produced in the B decay). For large values of 
this cut gets closer to another cut of the amplitude, corresponding to different physical processes, 
defined for v ■ q > {(2Mb + -^/i)^ ~ q^ ~ M^)/{2Mb). More precisely when takes its maximum 
value {Mb — Mh)'^, the first cut is reduced to one point v ■ q = {Mb — Mh) and the other starts 
aX, V ■ q = Mb + 3M/j. In our case M^ is the pion mass and we must expect the breakdown of 
the expansion when q^ is maximum. We remark that this problem is not related to the 5 and 5' 
in ( p!7| ) and ([l9|), those appear as soon as we allow k 7^ independently of the value of q^ . In 
our calculation no hadronic mass M^ appears and since all the decay products are massless, the 
maximum value of q^ is m? . In the electron spectrum we have q^ = m? at the endpoint Ee = E^"'^. 
Conversely, in the case of dTB/dw, we find q^ = rm? at the opposite edge of the spectrum that is 
at w = 0. 

From (|l7|) and (jl^), since is negative, we see that the spectrum takes negative values around 
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u; = in a region of the same order of magnitude as Gh- Although this is a very smah effect because 
of the value of Gb, it is clear that our spectrum is a poor approximation to the physical one which, 
as a matter of fact, must be zero for any value of the energy below the first threshold. 

With no information about the subsequent terms in the series, one can reasonably take h.QCD 
as an estimate of the size of the energy range where the expansion breaks down, and eventually 
consider the spectrum in (|T7|), (|l9|) as a realistic approximation for energies above, say, .8 Gev. 
Therefore, we suggest the possibility of using the hadronic spectrum to determine Vub- In fact b ^ c 
contributes to the B decay only for energies greater than the charm mass, so that b ^ u can be 
safely considered the only allowed channel up to approximately 1.1 GeV. This means that there is 
an energy interval left, of about the same size as the one available at the endpoint of the electron 
spectrum, where the hadronic energy distribution can be of interest, being free from the problem 
of singularities that affects the electron spectrum at Ef> = E^""^ . 

In conclusion we have evaluated the B decay spectrum as a function of the massless final 
quark energy at the 1/m^ level and to lowest order in Os/tt and we have found that it shows, at 
the maximum energy endpoint, the same kind of singularities as the electron spectrum and, at the 
same time, it turns out to be unphysical also in the low energy region that is typically dominated by 
non perturbative effects. Still it seems that there is a small range, usable for the Vub determination, 
where the approximation holds. 
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Figure Caption. 

Differential widtli dT / dw as a function of w as given by the sum of (|l^ and (H). The value of the 
parameters is Gb = —0.0065 and Kb = 0.01 (solid line), Kb = 0.0 (upper dotted curve), Kb = 0.05 
(lower dotted curve). 
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